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Abstract 

We investigate the phase transition in the three-dimensional abelian Higgs model for N 
complex scalar fields, using the gauge-invariant average action IV The dependence of 
Tk on the effective infra-red cut-off k is described by a non-perturbative flow equation. 
The transition turns out to be first- or second-order, depending on the ratio between 
scalar and gauge coupling. We look at the fixed points of the theory for various ./V 
and compute the critical exponents of the model. Comparison with results from the 
e-expansion shows a rather poor convergence for e = 1 even for large N. This is in 
contrast to the surprisingly good results of the e-expansion for pure scalar theories. 
Our results suggest the existence of a parameter range with a second-order transition 
for all N, including the case of the superconductor phase transition for N = 1. 
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1 Introduction 



The quantitative details of the phase transition in superconductors are so far not 
fully understood. If the phase transition is second-order or very weakly first- 
order, it is believed that the universal behaviour in the immediate vicinity of the 
critical temperature can be described by a field theory of electrodynamics with 
a charged scalar field. Here the scalar field may represent composite degrees of 
freedom like Cooper-pairs in "standard" superconductors. A non-vanishing ex- 
pectation value of this scalar field indicates spontaneous symmetry breaking and 
corresponds to the superconducting phase. Even though the critical behaviour 
for low temperature superconductors may be difficult to observe, this situation 
may improve for high temperature superconductors [I]. Also, it is claimed that 
certain phase transitions in liquid crystals are described by the same universality 
class as for the superconductors J|. Critical exponents have been measured for 
these transitions ||. An understanding of the critical behaviour in the three- 
dimensional abelian Higgs model, which should describe all these transitions, is 
also interesting from the theoretical side: This is the simplest example with a 
continuous gauge symmetry. Furthermore, the high temperature phase transi- 
tion in the four- dimensional abelian Higgs model should be well approximated 
in the vicinity of the phase transition by a three-dimensional effective theory [Q . 
Here the abelian gauge theory may serve as a prototype for the understanding 
of the non-abelian model which describes the electroweak phase transition in the 
early universe. 

In this paper we address the problem of the phase transition in the three-dimensio- 
nal field theory of scalar electrodynamics. Even the order of the phase transition 
is not settled so far, ranging from proposals of a first-order transition for all 
values of the coupling constants to speculations about a possible second-order 
behaviour |J . On the one side high temperature perturbation theory in the four- 
dimensional model seems to predict a first-order transition. Perturbation theory 
is, however, plagued by severe infra-red divergences. In the symmetric phase 
the one-loop correction to the quartic scalar coupling diverges as k" 1 , where 
k is an appropriate infra-red cut-off. Although it is possible to improve the 
calculations by a resummation of diagrams, there are questions on the reliability 
of this method [0]. Moreover, from the one- loop correction to the gauge coupling, 
an infra-red divergence ~ k^ 1 appears if the scalar mass vanishes. Since this 
is the case approximately near the phase transition, one concludes that in this 
region the results of perturbation theory will not be accurate. The encountered 
infra-red problems are exactly those of three-dimensional perturbation theory. In 
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order to avoid this difficulty an alternative approach uses the e-expansion and 
attempts to extrapolate results of a computation in 4 — e dimensions to e = 1. 
This has not given a convincing picture so far. The situation improves somewhat 
if, in addition, one considers N charged scalar fields with N large. The phase 
transition is well understood for N — > oo but an extrapolation to iV = 1 is not 
easy. We will see that even the leading 1/N contribution to the critical exponents 
is not well described by the e-expansion. 

In recent years, a new method which deals with these infra-red problems in a 
non-perturbative way has been proposed ||. The method is based on the ef- 
fective average action, r&, for which only fluctuations with momenta q 2 > k 2 
are included. The scale k acts as an effective infra-red cut-off. An exact, non- 
perturbative evolution equation describes the scale dependence of the average 
action [|]. Using this equation, it is possible to follow the /c-dependence of T k 
to extrapolate to the limit k — > 0, where Tk becomes the generating functional 
for the 1PI Green functions. The method has been tested in the $ 4 -theory fllO 



where it successfully describes the second-order phase transition that the system 
undergoes. A study of the critical behaviour of the three-dimensional $ 4 -theory 
gives values for the critical exponents in very good agreement with previous es- 



timates JTTJ. A gauge- invariant effective average action can be constructed for 



abelian and non-abelian gauge theories \12\. The corresponding non-perturbative 
exact evolution equation for the abelian Higgs model has been derived for arbi- 
trary dimension d For d — 4 and small gauge coupling, the flow equation 
is similar to the perturbative renormalization group equations in the one-loop 
approximation. A numerical study, based on this work, indicates a first-order 
character for the Coleman- Weinberg phase transition jn|, in agreement with the 
perturbative analysis. However, for d = 3, strong renormalisation effects of the 
gauge coupling result in predictions that may not be obtained by perturbation 
theory. In particular, the dimensionless coupling e 2 reaches in the scaling limit a 
value of order 6n 2 , thus e 2 /167r 2 is not a small quantity anymore i. 

The order of the phase transition seems to depend on the ratio Ao/e 2 , of the 
bare quartic scalar coupling Ao and the squared gauge coupling e^. For small 
Xo/ e o (tyP e I superconductors) it was suggested that the phase transition is first- 
order H, the discontinuity being induced by fluctuations of the gauge field. This 
picture has been confirmed by the solution of the non-perturbative flow equation 
HJ. The question about what happens for large scalar couplings has not found 



5 For the non-abelian Higgs model in three dimensions the dimensionless gauge coupling also 
grows large. Additional complications arise from confinement effects 03. 
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a clear answer by previous methods, except for very large N (typically iV > 200) 
where a second-order phase transition was established for sufficiently large A /eQ. 
Investigations based on the average action indicate the existence of a range in 
parameter space where the transition is second-order also for the superconductor 



(N = 1) If a parameter range for a second-order transition exists for all N 
there must be a fixed point of the flow equations with associated scale invariant 
critical behaviour. 

In this paper we will use the concept of the average action to examine the critical 
behaviour of the three-dimensional abelian Higgs model for N complex scalar 
fields. We will discuss how the fixed points of the theory move as iV varies and 
examine the critical exponents of the model, in the region where a second-order 
phase transition arises. In section 2 we give a brief sketch of the basic formalism 
and derive the evolution equations for the running of the couplings. In section 
3 we look at the scaling solutions in the large-iV approximation. In section 4 
we discuss a differential equation which allows us to follow the iV-dependence of 
the fixed points. This uses a linearization of the flow equation for values of the 
couplings near the fixed point values. In section 5 we give a detailed numerical 
analysis of the model for various N. The critical exponents for a second-order 
phase transition are discussed in section 6. In section 7 we present our results and 
compare them with those obtained through the e-expansion. Finally in section 8 
we give a summary of the work and derive the conclusions. 



2 Average action and flow equations 
for running couplings 

The starting point of the analysis is the classical action S for N complex scalar 
fields Xa, a = 1, ..,N and the gauge field A^ 

S[ X ,A,] = (1) 

where = d^A u — d u A^, = <9 M + ieA^. We note that the pure scalar 
sector is 0(2iV)-symmetric. The coupling to the photon reduces this symmetry 
to SU(N) x U(l) where the abelian part is gauged. Furthermore, we have the 
gauge fixing term 

S gf = ~7^ J d d x(A^ - A^d v {A v ~ A"), (2) 
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where A^ is some arbitrary background field. The average action Tk is the effective 
action for (gauge-invariant) averages of fields. It is obtained by integrating out 
only quantum fluctuations with (covariant) momenta larger than k. For this 
purpose one introduces in the functional integral for the usual effective action 
an infra-red cut-off Rk- If Rk vanishes for k — > one recovers the generating 
functional for the 1PI Green functions T as lim^o IV The functional Yk[tp, A] 
is gauge-invariant for all k. The dependence of the effective average action on 



the scale k is determined by an exact non-perturbative flow equation [jrj] which 
reads (up to minor modifications) 

|r t [„ , A] = ba (r»[„ , Al+vfj+R^y^lA]. (3) 

Here the trace involves a momentum integration and summation over internal 

(2) 

indices, t = Ink/ A, and Y k is the second functional variation of IV The sum 

IV + Y j corresponds to to the exact inverse propagator in presence of the fields 

ip and A. We approximate Y g f by (§) with replaced by A^. More details can 

be found in ]I2] where also the correction Cfc[A] is motivated. 

The infra-red cut-off Rk acts on scalars (R k ) and gauge fields (R k ) and reads 

explicitly 

Rl S) = Z,, k .(-D\A)^ ti^ D2 ^ 



fli-D\A)) 
fl i-d 2 " 



Here and Zp k are appropriate wave function renormalization constants for 
the scalar and gauge field, respectively. We choose a cut-off that decays expo- 
nentially for large (covariant) momentum squared x = q 2 , i.e. 



fl{x) = e-'". (5) 



For x — > the cut-off R k — > Zk 2 becomes a mass-like term and suppresses the 
propagation of quantum fluctuations with q 2 < k 2 . This eliminates the infra- 
red-divergence problem that appears in the perturbative approach. For large 
q 2 the exponential decay of d t R k guarantees the ultra-violet convergence of the 
momentum integral in (§). Finally, the use of the covariant derivative in (f|) 
ensures explicit gauge covariance of the flow equation. 

A discussion of the phase transition relies mainly on the properties of the effective 
scalar potential U k which corresponds to the free energy. In order to extract from 
(H) the flow equation for the average potential Uk we use a configuration with 
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ip = const, and A = 0. For the derivation of the flow equation for U k we use a 
truncation of the average action of the form 

r k [<p, A} = J d d x (^Z Fjk F^ + Z v>k (D^ a )*D^ a + U k {p)) . (6) 

The second functional derivative of this expression can be found in appendix A. 
One obtains for arbitrary dimension d the flow equation 





Q-Uk(p) = (d- l)v d J dxxi-% In (Z F>k P(x) + 2Z^ k e 2 p 

+v d J dxxt-% In (Z^ k P(x) + U' k (p ) + 2p U' k \p )) 

+ {2N - l)v d [ dxxi-% In (Z V)k P(x) + U' k {p )) , (7) 



where 

w -i = 2 *V^). (8) 

Here the potential is a function of the invariant p = f a f a and primes denote 
derivatives with respect to p. We also use the variable x = q 2 and the inverse 
average propagator P(x) is defined by 

ZP{x) = Zx + R k (x) 
1 - exp(- F ) 

The partial derivative dt acts only on the infra-red cut-off R k in P with 

d t P(x) = k^P(x) - (k^-lnZ k )(P(x)-x) (10) 

and we drop the last term proportional to the anomalous dimension in this paper. 
In the sequel, the use of dimensionless variables turns out to be convenient. We 
therefore introduce the following variables: 

u k (p) = k~ d U k (p) (11) 
p = k 2 - d Z^ kP (12) 
e 2 = fc d - 4 (Z J , fc )- 1 e 2 . (13) 

In terms of the above, the partial differential equation (|7]) can be rewritten as 

d t u k (p) = -du k (p) + (d-2 + r] v )pu' k (p) + 2(d-l)v d l d s d (2e 2 p) 

+2(2N - l)v d l d s d (u' k (p)) + 2v d l d s d (u' k (p) + 2p<(p)), (14) 
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where primes denote now partial derivatives with respect to p. The anomalous 
dimension r) v = —d t In reads 

= -16v d (l - ^) e 2 /f,i(2A/t, 2e 2 n) + ^jv d \ 2 n mJ 2 (2A/t, 0) , (15) 

and we introduced the dimensionless renormalized vacuum expectation value 
k defined through u' k (K) = and the dimensionless renormalized quartic cou- 
pling A as 



X(k) = <(«(*)) . (16) 

The constants l d and the threshold functions Sq(lj), l d m (uJi, u 2 ) and m d m (uJi, uj 2 ) 
are given in appendix B. We also need an evolution equation for e 2 (k) that we 
infer from Jl2| . 



The main subject of the present paper is the investigation of solutions to eq. 
In order to solve eq. (|14]) for Uk{p) we make first an Ansatz for the effective 
potential Uk(p) as a polynomial expansion around the /c-dependent minimum 
at K,(k). The partial differential equation ( |14|) can then be transformed into 
infinitely many coupled ordinary differential equations for the higher derivatives 
of the effective potential taken at the vacuum expectation value. To lowest order 
we find the following coupled set of differential equations for the regime with 
spontaneous symmetry breaking 

de 2 

" & 2 = (d-A)e 2 
4 



dt 



+ rd e* 



lp d g (2\K, 2e 2 n) + l d c (2\K) + (N- 1)1% (17) 



^= fi K = (2-d- Vip )K + 4j(d-l)v d l d s d (2e 2 K) 



d\ 
~dt 



+6v d l d s d {2\K) + 2{2N - l)v d l d (18) 
Px = (d - 4 + 2r^) A + 8e 4 (d - 1) v d l d s d 2 (2e 2 K) 

+18\ 2 v d l d s d 2 {2\K) + 2(2A^ - 1)A 2 v d l d . (19) 



The functions s d (u), l d {uj) and Sg(ui,Uz) describe the threshold effect due to 
mass terms. They are given, as well as the constants l d , l d c and l d , in appendix B. 
We emphasize that our /^-functions are directly computed for arbitrary dimension 
d and can immediately be evaluated for d = 3. No e-expansion around d = 4 
is necessary. We will see later (sect. 7) that the e-expansion for e = 1 fails to 
reproduce our results. 
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The iV-dependence in f3 K and f3\ is coming from diagrams that involve the (2N— 1) 
massless scalar modes inside the loop. We note that the scalar contributions 
are the same as for a pure SO(2N) -symmetric scalar theory for which one has 
(2N — 1) massless Goldstone excitations around a minimum at k ^ 0. For f3 e 2 
the iV-dependence arises from the diagrams that involve the (N — 1) fields a 
and uj (appendix A). They correspond to iV — 1 complex massless scalars which 
have a vanishing expectation value. These 2N — 2 real Goldstone modes be- 
long to the fundamental representations of the unbroken global symmetry group 
SU(N — 1). The contribution of the remaining "would be" Goldstone boson 
(singlet of SU(N — 1)) is contained in the term for N — 1. 

Starting with sufficiently large initial values of A/e 2 and solving the system of 
equations (|17|)-(|T9|) for k — > shows that k either runs to zero at some non- 
vanishing scale k s or diverges ~ fc -1 . The first behaviour corresponds to the 
symmetric phase where po(k = 0) = whereas the second indicates spontaneous 
symmetry breaking with p (k = 0) = po > 0, k = pok' 1 . By an appropriate 
tuning of the initial value of k the scales k s or po can be made arbitrarily small. 
The phase transition is second-order. For the initial n taking exactly the critical 
value, the solution of the flow equations runs asymptotically towards a scaling 
solution for k — > 0. This is described by fixed point values k*, A* and e 2 and 
corresponds to the scale invariant physics at the phase transition. On the other 
side, for a small ratio of initial values A/e 2 the coupling A may reach zero (or 
extremely small values) for nonvanishing k. Looking at the flow equation fll4D 
for Uk this corresponds to a situation where already the absolute minimum of 
has jumped to the origin at p = 0. This behaviour indicates a first-order phase 
transition. 



3 Scaling solutions and large- N approximation 

Thus, we find that the abelian Higgs model in three dimensions has a region of 
parameter values where the phase transition is first-order, and a region where 
it is second-order. This can be qualitatively understood by looking at the two 
limits of the theory: When no gauge-field fluctuations are present (e 2 = 0) the 
model reduces to the 2iV-component Heisenberg model with a second-order phase 
transition. Below the critical temperature the effective potential Uq has a max- 
imum at p = whereas the minimum occurs for p = Po{k — 0). On the other 
hand, gauge- field fluctuations dominate the evolution equations for A -C e 2 . The 



solution of the flow equation is then similar to the four- dimensional model [13 
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and results in a first-order phase transition. A local minimum remains at the 
origin even somewhat below the critical temperature. 

At the critical temperature T c of a second-order phase transition the theory be- 
comes scale invariant. This corresponds to a fixed point where the dimensionless 
parameters k, A and e 2 take fc-independent values. In consequence, the unrenor- 
malized quantities po, U^(po) and e 2 vanish with appropriate powers of k. This 
"second-order fixed point" is infra-red stable in all couplings except one "relevant 
coupling" that we may associate with k. We specify short distance (or micro- 
scopic) couplings at some appropriate high momentum scale A. Then the critical 
surface, for which the phase transition occurs, is given by a critical value k c (A) 
which is a function of A(A) and e 2 (A). The difference k(A) — k c (A) is propor- 
tional to T c — T. With initial values on the critical surface all couplings flow 
towards the second-order fixed point if one starts within its range of attraction. 
This range of attraction also defines the region in parameter space for which the 
transition is second-order. There must also be a line (more generally a hypersur- 
face) which separates this parameter region from the one for which the transition 
becomes first-order. On this line, the couplings flow towards another fixed point, 
the "tricritical fixed point". This tricritical fixed point has two infra-red unsta- 
ble directions (in k and A). Our picture of a parameter range with second-order 
transition and a different range with first-order transition is thus realized if two 
non-trivial fixed points can be identified, one with one and the other with two 
unstable directions. 

We therefore want to investigate the scaling solutions related to the fixed points 
of the coupled system of differential equations (|TT|) - (O) in three dimensions, i.e. 



we look for simultaneous zeros of the r.h.s. of eqs. (|PTD - (pi). It turns out that 



the system ([17]) - (|19|) has non-trivial fixed points for any N, although previous 
results 0, obtained through an e-expansion method, seem to indicate that a 
large number of complex fields would be needed for this purpose. The main new 
ingredient for the existence of fixed point solutions even for small N stems from 
the fact that threshold effects due to the decoupling of massive fluctuations are 
properly taken into account. As an illustration i , we have depicted in fig. 1 the 
/5-function for A fll9|) as a function of A at fixed e 2 = e 2 . Neglecting the threshold 
effects corresponds formally to k = and results in the upper curve. In this 
case, no fixed point for A is obtained. Taking the full threshold behaviour into 
account (with k at its fixed point k+) reveals the occurrence of two fixed points 



6 We have chosen N = 10, but fig. 1 describes the generic behaviour for any N below ~ 183. 
For N > 183, fixed points are present even when the threshold effects are discarded. 
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for A (zeroes for (3\). Here the smaller value of A corresponds to the tricritical 
and the larger one to the second-order fixed point, respectively. 

To start a quantitative discussion of the location of the fixed points we consider 
in this section the large- N limit. The evolution equation for the gauge coupling 
decouples now completely from the scalar sector and reads 

Its non-trivial infra-red fixed point solution is given by 

e * = (21) 

The corresponding solutions for the scalar sector can be obtained through an 
Ansatz for the large- N behaviour of the couplings as k ~ N a , A ~ N 13 . Inserting 



this Ansatz in eqs. (|l8f ) and (|T9|) it follows that only the pairs 



(«,/?) = (1,-1); (1,-2) (22) 

may yield solutions to the fixed point equation. This result obtains also through 
an appropriate rescaling of the full partial differential equation ([14]) and is there- 
fore no mere effect of the polynomial approximation. Here A* ~ 1/N corresponds 
to the second-order fixed point, whereas the tricritical point is characterized by 
A* ~ 1/N 2 . More explicitely we find for the second-order fixed point (j3 = — 1) 
that the influence of the gauge coupling (^l|) and of the massive scalar fluctuation 
in the evolution equations (|1|) - ([19]) are of sub leading order in 1/N. Therefore, 
they simplify to 

_ = - K + 4 Nvsl * (23) 

§ = -A + 4iV^A 2 , (24) 
at 

with the corresponding fixed point values given by 

k* = Avzl\N (25) 

= 555 F (26) 

Since eqs. (|2H|), (^3]) and ( PH ) are decoupled, it is obvious that the fixed point 
characterized by eqs. (|2l]) , (|25|) and (|26|) is infra-red stable in A and e 2 and unstable 
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in k. Since it has only one unstable direction this must be the second-order 
fixed point. The fixed point values k+ and A* are entirely determined through 
the fluctuations of the massless modes and any memory of the presence of the 
gauge field has disappeared in this limit. Therefore, in the large- N limit, the 
fixed point values A* and are identical to the Wilson-Fisher fixed point of the 
0(2iV)-symmetric pure scalar model. This simple relation does not hold true 
for the subleading terms in A* nor for the critical indices of the theory. The 
anomalous dimension, for example, reads for the fixed point (pT|) , fl25|) and (|26|): 



4 If 



m 



3 f™2 
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b 2) 



v 



i 



(27) 



where m* (resp. M*) denotes the dimensionless scalar- (gauge-) field mass at the 
fixed point: 



ml 



Ml 



2/ 3 

2A*re* = = 1.7 
h 

2eX = -i = 6.3 . 

l gc 



(28) 
(29) 



The first term in eq. (27) stems from the scalar fluctuations only (and corresponds 
to the known result for the pure scalar case), whereas the second term describes 
the effect of the gauge-field fluctuations. Indeed, the latter one is the dominating 
term and we find that rj* eq. (p7]) is negative 



V* 



0.31 



(30) 



Next we turn to the second solution given by (3 = —2. Here, the situation changes 
significantly since eq. (|i~9l) simplifies to 



~dt 



-A + 16e% 3 llsl(2e 2 K) 



(31) 



The fixed point is now given by 



A* = 16e^ 3 ^4(2e;fre 



9/3 2 



1 



(32) 



In this region, the leading contributions to eq. Q3~T| ) are coming solely from the 
gauge-field fluctuations (~ 1/A^ 2 ), and even the massless scalar fluctuations con- 
tribute only in subleading order (~ 1/iV 3 ). The fixed point ( |3"2"D is thus a genuine 
effect of the presence of the gauge field, not present in the pure scalar case. For k± 



10 



we find a fixed point ~ N which corresponds to the solution of the A-independent 
implicit equation for 2e 2 A* 



3^2 



2e*K* = 8v 3 qe 



A + 2^(2e^ 



6g t{ s\{2el^) 



Comparing eq. (|29f) with eq. ([33|) we see that the numerical value for n+ is affected 
by corrections from the gauge-field fluctuations. (Remember that the fixed point 
e 2 . eq. ([H]) is the same for both solutions eq. (p2p.) We wish to put emphasis on 
the fact that this fixed point is strictly different from the Gaussian fixed point 
even in the limit iV — ► oo. Analyzing the stability properties of the fixed point 
specified by (j21[) , ( 32|) and (|33|) reveals two unstable directions (k and A). We 
associate this point with the tricritical fixed point. 

In conclusion, we have found in the large- A approximation two non-trivial fixed 
point solutions of eqs. (jl7|)-(|T9l). One solution governs the second-order phase 
transition in this model. The second fixed point is mainly an effect due to the 
gauge-field fluctuations and governs the tricritical behaviour. Starting, for exam- 
ple, with ka near re c (A), e\ = e 2 (A) but A (A) < A* (A) (cf. eq. (|32|) ) the model 
undergoes a first-order phase transition as k(A) passes through k c (A). Additional 
fixed points occur for e 2 = 0: There is the infra-red unstable Gaussian fixed point 
at A* = with k+ given by eq. fl2"5D and the Wilson-Fisher fixed point of the pure 
scalar theory with A* and given by (|26D and (f2~5l). Since e 2 = this fixed point 



is different from the second-order fixed point and has two unstable directions. 
Finally, if we consider 1/e 4 instead of e 2 two more similar fixed points exist for 
(l/e 4 )* = 0. Knowledge of all the fixed points and their stability properties per- 
mits easily to establish the qualitative features of the whole phase diagram. 
In section 6 we will give a detailed numerical analysis of the second-order fixed 
point and the tricritical fixed point for arbitrary N. In the same section we will 
also extend the analysis beyond the very simple approximation (fl7D - (fL9|) and 
discuss more elaborate approximations to the flow equation eq. (]14]) . 



4 TV-dependence of the fixed points 

For large N we may obtain the fixed points of the theory by just looking at the 
leading A-dependent contributions in the beta functions. The A-dependence of 
the location of the fixed points can be discussed in terms of a simple differential 
equation which involves the derivatives of the /3-functions at the fixed point. A 
solution of this differential equation, with initial conditions set at large N, allows 
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to compute the fixed points also for small values of N. In order to derive the 
flow equation for the iV-dependence of the fixed points, we define the generalized 
vectors 

/ k \ ( d t K \ 



A 



A 



d t \(N)={3(\(N),N) 



(34) 



(This discussion can be easily extended if additional couplings are taken into 
account). In general, A is a function of k and N. However, since we are interested 
in the scaling solutions for arbitrary N, we wish to take A on a fixed point. 
Therefore, every scale dependence is removed and A becomes solely a function of 
N, \+(N). For all iV it obeys the fixed point condition 



dt\*(N)=/3(K(N),N) = 0. 



(35) 



Taking the total derivative of eq. (|35|) with respect to iV gives the desired differ- 
ential equation for the iV-dependence of A*: 

^(N) = -A- 1 (A,(iV),iV)|^(A,(iV),iV). (36) 
Here, the matrix A(X(N), N) is the Hessian matrix at the fixed point, 

A(\,N) = ^(\,N). (37) 

CM 

For given N, it governs the evolution of linearized fluctuations around the scaling 
solution 



dt(\-K)=A(\-K) ■ 



(3* 



As long as A remains regular, the flow equation has a solution and the fixed 
point exists for all A" in this domain. Also the stability property of the fixed point 
cannot change in the domain where A remains regular. We have plotted in fig. 2 
the eigenvalues of A for the second-order fixed point in dependence on N. All 
eigenvalues remain far away from zero and depend only moderately on N. Within 
our truncation we conclude that a second-order fixed point and the associated 
parameter region for a second-order phase transition exists for arbitrary N > 1. 
(This can formally also be continued to AT < 1). We have numerically solved the 
flow equation (|3B|), starting initially with large A^ and extrapolating to A^ = 1. 
The results that we obtained using this method are in very good agreement with 
those of the complete numerical analysis that we present in the following section. 
The difference between the two methods is indistinguishable in the plots. 
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5 Numerical Analysis 



In this section we present in detail the results of a numerical investigation of the 
fixed-point structure of the abelian Higgs model. Since we are especially inter- 
ested in the critical behaviour of the theory, we will focus on the existence of 
fixed point solutions and the related critical indices. As we have mentioned, we 
employ an expansion around the "asymmetric minimum" at k ^ 0, thereby trans- 
forming the partial differential equation in an infinite series of coupled ordinary 
differential equations for the corresponding couplings. This system, then, has to 
be truncated at some finite order and will be solved numerically. 

In lowest approximation, the differential equations are given by eqs. (|17D - (|19|). 
The numerical values at the second-order fixed point of n+ , A* and are given 
as functions of the number of complex scalar fields N. Our results are depicted 
by diamonds in figs. 3-5. The solid lines correspond to the large- N extrapolation 
as given by eqs. fl2"T|), and (PB|). It turns out that is very well described by 
eq. ( Pq) even for N = 1. This is related to the fact that the contributions from the 
massive fluctuations in eq. fll8|) are damped through the threshold behaviour. The 
subleading term at AT = 1 stems from the third term in eq. ( |18|) and is suppressed 
by a factor of approximately 3sf(2A*/c*) ~ O(10~ 2 ) compared to the leading one. 
The contribution of the gauge-field fluctuation in eq. (|T8| ) remains unimportant 
despite the rather large value for the fixed point of the gauge coupling (e| ~ 633 
for iV = 1). It is suppressed by a factor of 4sf (2e^K^)e^/A^ ~ C(10~ 3 ) compared 
to the leading term. Fig. 4 displays the results for A* . Again, we find that A* is 
very well described by its large- N extrapolation eq. (f26"|). Only for JV = 1 we 
observe a sizeable correction due to the third term in eq. (|19j) which becomes 
nearly of the same order as the last one. Their ratio reads 9sI(2\+k+) ~ (9(10 _1 ) 
and small corrections are of no surprise. The contribution of the gauge-field 
fluctuation in eq. ( ]T9| ) is suppressed by a factor of 8sl(2elK*)ey \l ~ O(10~ 3 ) 
and therefore again unimportant. The results for el are given in fig. 5. The 
large-iV extrapolation eq. is in a very good agreement with our numerical 
results down to iV ~ 10. For iV < 10, approximating all fluctuations as being 
massless becomes less and less accurate. The mass of the fluctuations will lower 
the e 4 -coefficient in eq. ( |TTD as compared to the purely massless estimate and 
therefore enhance the numerical value for el . (Clearly, this effect is largest for 
N = 1 where all contributions in eq. (|17D are suppressed by massterms and the 
estimate eq. (|2"T|) becomes too small by roughly a factor of 10.) 



In order to check the stability of these findings with respect to a change in the 
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truncation, we studied approximations of the effective potential being a local 
polynomial in (p — k) up to (p— k) 3 and (p — k) 4 . (We call these approximations 
the <p 6 - and <p 8 -approximation, respectively). This approach turned out to be 
very successful in the pure scalar model [10] , where the non-trivial scaling solution 



corresponds to the Wilson-Fisher fixed point. However, we do not expect that this 
approach will give reliable results near the tricritical fixed point. Quite generally, 
the local polynomial approximation around the asymmetric minimum at po(k) 
breaks down when X(k) reaches zero for k > in the course of its running. This 
happens near a first-order phase transition and is connected with the fast running 
of k and all the higher couplings in the region of very small A as can be seen from 
eq. (|I%D. A discontinuity in k (jumping from k c > to k = at A = 0) will be 
manifest through a pole ~ 1/A in the /^-functions. Nevertheless, these regions in 



parameter space can be handled by a more global approach as described in [13 
for the d = 4 abelian Higgs model. Here, flow equations at the local minimum 
and at the origin are used simultaneously. At the two extrema we expand in 
second-order in p or (p — p ) and guarantee continuity of Uk by assuming a p 4 
polynomial for the region between and po- By this method, the evolution of po 
knows about the flow at the origin (for example the appearance of a new local 
minimum) and vice versa. As long as Uk exhibits a local asymmetric minimum 
(A > 0,k > 0), the effective average potential is parametrized in terms of the 
four variables A, K,m^/fc 2 , and X s , where X s (k) = Z~ 2 k U' k '(p = 0)/c _1 denotes the 
dimensionless quartic scalar self-coupling, and m 2 s {k) = Z Vi jJJ' k {p = 0) the scalar 
mass term at the origin. This second method will serve as an alternative check 
of the results from the local polynomial approximation at the second-order fixed 
point in the various approximations. In addition, it will yield results for the 
tricritical fixed point. 

It turns out that is very stable against changes in the truncation. The results 
for the <p 6 - and <p 8 -approximation and those obtained with the global method 
are essentially identical to fig. 3 and differ only slightly from them for small N. 
The fixed points for e 2 and A* within the different approximations are depicted 
in fig. 6. It shows e 2 as a function of A* for different values of N. The solid line 
corresponds to the large- N estimate eq. (pTf) and eq. (|26|): 



3/3 

e'(A*) = -^A* = 3.7A* . (39) 

l 9C 



All the different approximations converge rather fast to the large- N result eq. fl39|) . 
However, for small values of N the fixed points are quite different. But this has to 
be expected: The inclusion of a coupling ~ <p 6 (being dimensionless in d = 3) is 
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for small iV an important effect for A* [11]. In turn, depends now strongly on 
A* since the massive excitations are important, which results in a sizable shift of 
the second-order fixed point. On the other side, we do not expect that couplings 
~ ip 10 will modify the result of the ^-approximation in an important way. In 
summary, the second-order fixed point is confirmed in a variety of approximations 
to the exact evolution equation. The estimates eqs. (PTJ), and (BB) give very 



good predictions already for small N. The results for the critical exponent v and 
the anomalous dimension rj and their discussion is given in section 7. 

Finally, we turn to the tri critical fixed point. In the ^-approximation we find a 
tricritical point for all N > 2. For N — 1 the gauge coupling runs to infinity in the 
interesting region and the tricritical point presumably corresponds to e\ — > oo. 
This may well be an artefact of our truncation. Our results obtained by the global 
approach and the (^-approximation are illustrated in fig. 7. Again, e\ is given as 
a function of A* for various values of N. The solid line represents the expected 
behaviour for large N as given by eq. (|2T]) and eq. fl32|): 



<(A*) ~ y/K- (40) 
This is well reproduced by the global method. Here the results from the <p 8 - 



approximation converge rather slow to the asymptotic behaviour (|40|). A com- 
ment on the solutions for small N is in order. For N < 4 no tricritical fixed point 
has been found with the global method. This is related to the fact that we have 
to expand the potential at p = around a local maximum. The running of the 
dimensionless mass term at the origin 

d t m 2 s = -2 m 2 s - 8v 3 lle 2 + 0{X S , rj^ml) (41) 

is dominated for small N by the e 2 -term. The fixed point for ml reads therefore 

(ml)* ~ -4^4 (42) 

Requiring (ml)* to stay away from the pole of the threshold functions at ml = —1 
yields with the help of eq. (pi]) the condition 

N > |i = 3.2 , (43) 

which explains roughly the absence of fixed points for N ^> 4 in this approach. 
However, we do not expect this to be a physical effect. Instead, it is related to the 
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fact that we used in the global approach the evolution equation QT7| ) not only at 
the asymmetric minimum, but also at the origin. This approximation is harmless 
in the large- N limit where el (eq. PIP) and the fixed point for the dimensionless 
renormalized gauge coupling at the origin e|j* are both small. Their difference 
behaves like 1/N. For small N, in contrast, this may result in a sizeable effect. 
As a result, is always smaller than el and the bound eq. may therefore 
be lowered easily.1 

In conclusion, it has been shown that the second-order and tricritical fixed points 
are confirmed numerically in a variety of different approximations, locally and 



globally, to the exact evolution equation fll4l). The second-order fixed point is 
shown to be rather stable against higher truncations, especially for large N. For 
large N we expect a rather fast converge of the series of approximations <£> 4 , 
ip 6 , ip 8 ... towards the scaling solution of the partial differential equation (|14]) as 
given by d t Uk(p) = 0. Since all contributions of terms neglected in the truncation 
(^) are suppressed by powers of 1/N, the predictions of critical exponents com- 
puted from this fixed point should be very reliable. For small N the convergence 
of the series of truncations is not obvious with the restricted Ansatz ([]). The 
reason leading to the "missing tricritical point" in the global method will also 
affect the existence of a scaling solution of eq. (|14]) and therefore the convergence 
of the series. For small iV it is not excluded that the (^-approximation will give 
better critical exponents than higher truncations. For the tricritical point the 
region of convergent truncations is restricted to even higher values of N. 



6 Critical exponents 

The behaviour of a system near the critical temperature of a second-order phase 
transition is described by the critical exponents of the theory [TjJ . These specify 
the long range correlations of physical quantities as we approach a zero mass 
theory at T = T c . The system is described by five exponents, which however are 
not independent. For example the indices 7 and S, which describe the response of 
the vacuum expectation value p to an external magnetic field or source, as well as 
(3, which specifies how the vacuum expectation value depends on the temperature, 
may be calculated, once both temperature dependence of the correlation length 
and the behaviour of the connected two-point function at the critical temperature 
are known. Thus, we will look how the two independent indices v and 77 behave 



7 In order to make a full account of this effect, the p-dependence of Zp(p) (replacing Zp in 
eq. (0)) has to be calculated. We refer this investigation to future work. 
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as N varies. How these can be computed from properties of the scaling solution 
is the subject of this section. 

The critical exponent r\ determines the behaviour of the connected two-point 
function at the critical temperature (T = T c ). 

lim < <p(x)(p(0) > c ~ \x\- {d - 2+Jl) . (44) 

x — >oo 

It arises from the anomalous dimension of the field ip and is directly related to 
the momentum dependence of the (unrenormalized) Green function for q 2 — ► 

at T = T c 



G-\q)=q 2 Z{ q ) (45) 



lim Z(q) ~ (q 2 )-. (46) 



The Green function G~ 1 (q) is encoded in the momentum dependence of the 
derivative terms in 1^ which are quadratic in small fluctuations around the 
if = const configuration. Its computation goes beyond the truncation (|6|) and 
would necessitate a generalization from a constant to a momentum depen- 
dent function Z k (q). Then Z(q) in eq. (flop obtains as lim fc ^ -Zfc(g) whereas Z Vs k is 
given by lim g 2^ Zk(q). Even without computing Zk(q) explicitly, simple scaling 
considerations indicate for the critical behaviour 




Here A is some appropriate high momentum scale and the smooth function / 
approaches constant values for both limits q 2 <C k 2 and q 2 ^> k 2 . This follows 
from the observation that q 2 and k 2 act as independent infra-red cut-offs in the 
loop describing the exact flow equation (|3|). In consequence, the critical exponent 
1] can be extracted from the fc-dependence of the wave function renormalisation 
for the scaling solution 



(4£ 




This allows to identify rj with rj v (15), where all couplings are taken at their fixed 
point values. 
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We now turn to the critical index v, which is related to the temperature depen- 
dence of the correlation length near the critical temperature T c 



lim < <p(x)(p(0) > c ~ e~ x/ ^ £ = m' 1 ~ |T - T c \~ v . (49) 



Here the renormalized mass m 2 is given by the limit k — > of the fc-dependent 
mass term m 2 {k) = Z^ k (k)U' k (0) and 2Z^ k U' k \p {k))p {k) = 2\{k) K {k)k 2 , for 
the symmetric phase and the phase with spontaneous symmetry breaking (SSB), 
respectively. The phase transition corresponds to a critical trajectory within the 
SSB regime ml(k) = 2\±K,±k 2 . 

In order to extract v we study a small deviation from the critical trajectory 

m 2 (k) = ml(k) + 5m 2 (k). (50) 
The evolution equation for 5m 2 is characterized by the anomalous mass dimension 

U! 

d 

—5m 2 = uj 8m 2 . (51) 
dt v ; 

For very small Sm 2 /k 2 we linearize in 5m 2 so that to becomes independent of Sm 2 

dm 2 . f3\ @ K 2 
/3 m 2 = = (2 + — + — )m . (53) 
at A k 

At the fixed point the anomalous mass dimension is a constant to* = u}+(k+, A*, e 2 ). 
We therefore obtain in the vicinity of the scaling solution for 5m 2 <C k 2 

5m 2 {k) = (j^j 5m 2 {k 1 ). (54) 

(Here ki denotes some scale below the ultra-violet cut-off for which all couplings 
are near their fixed points. It can be taken proportional to T c ). For uj < 2 the 
ratio 5m 2 / k 2 increases as k becomes smaller. For any non- vanishing 5m 2 (ki) there 
is necessarily a critical scale k c where the linearization of (|51~1) becomes invalid, 
typically for some constant c < 1 

bm\ = 5m 2 (k c ) = ck 2 c . (55) 
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Near the phase transition the renormalized mass is the only scale present. In the 
SSB phase we conclude from dimensional analysis that it must be proportional 

to k c 



m 2 = m 2 (0) = <W (0) = ak\ 



(56) 



The situation is similar in the symmetric phase, but the proportionality constant 
a is in general different for the symmetric and the SSB phase. We finally observe 
that m 2 (ki) is a function of temperature and that in linear approximation 



m 2 {h, T c ) = ml{h), Sm 2 (k 1 ) = b(T c - T). 



(57) 



(Remember that the symmetric phase corresponds to m 2 (ki) < m 2 (ki)). Com- 
bining eqs. (p4[),(|55D and (|57| ) gives 



and we conclude 



m = a 



\T-T e \r- 



(59) 



This is to be compared with eq. (f49|), showing that the critical exponent v is 
directly related to u+ 



1 



v 



(60) 



We emphasize that (0) with flBCf ) becomes exact in the limit T — > T c since the 
fixed point behaviour completely determines the evolution of m 2 . In other words, 
the "initial" and "final" running of m 2 away from the fixed point becomes neg- 
ligible compared to the running near the fixed point. The preceding discussion 
has been explicitely verified by numerical studies for the pure scalar theory [[IT]]. 
We still have to specify how the anomalous mass dimension should be evalu- 
ated. In fact the ratios ^ and ^ in eq. ( |53D depend not only on \k or m 2 which 
we have up to now taken as a variable for simplicity of the presentation. They also 
involve A and e 2 (and in principle also some other parameters characterizing the 
average action, like for example the <£> 6 coupling). For a computation of oj* we 
have to specify how the couplings we are using depend on each other. Instead of 
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m 2 we use here k as the independent parameter and express A and e 2 as a func- 
tion of k. (This can easily be extended if additional couplings are considered). 
Then u± can be expressed in terms of the partial derivatives of /3 K evaluated at 
the fixed point: 



2 + 



8k 



+ 



d\ (dp. 
Ok \ OX 



+ 



dn \ de 2 



(61) 



9Hk) QT1 ^ de (k) rj^ e p^gg transition 
corresponds to a running along a trajectory on the critical surface into 



It remains a determination of the quantities and 



(T = T C 

the second-order fixed point. Trajectories infinitesimally close to the critical 
surface (for T infinitesimally close to T c ) come infinitesimally close to the fixed 
point and then move away from it infinitesimally close to a curve denoted as 
"unstable direction". The dependence X(k) and e 2 (/t) on the curve corresponding 
to the unstable direction should be used in the equation for u;*. More formally, 
the evolution equation for small deviations from the fixed point 



5k = k — k+, 5X = X — A*, 5e 2 = e 



2 e 2 



(62) 



can be written as 



d_ 

dt 



( Sk \ 




( Sk\ 


5X 


= A 


5X 


V 6e2 ) 




\ Se2 ) 



(63) 



with (cf. sect. 4) 



A 



( dpjdn dpjdx d(3 K /de 2 \ 

df3 x /dK d(3 x /dX d(3 x /de 2 
K d(3 e 2/dK d(3 e 2/dX d[3 e2 /de 2 j 



(64) 



At the second-order fixed point all eigenvalues of the matrix A, except one, are 
positive. Let us denote by x us = (1, 62,^3) the eigenvector corresponding to 
the negative eigenvalue. As k decreases, only this eigenvector plays a role. It 
indicates the unstable direction whereas the components corresponding to positive 
eigenvalues quickly die out. For the unstable direction we therefore find that 
^ = 62 ; or, for the equation for u;*, |^ = 62, ^7 — &3- Since in general A is 
not a symmetric matrix, we should specify what we mean by an eigenvector. 
We denote by 9 the unique, negative eigenvalue of A. Using the transformation 
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A = C 1 AC we can bring A to a form where An = 65n. Then the unstable 
direction corresponds to the vector 

= crfxf = cS a (65) 

with 

^-x us = Qx us . (66) 

dt v ; 

Inverting (|65| ) we can read off b{ = where the coefficients Cn obey A^Cj\ = 
QCi\. This completely specifies to* by the derivatives of the /^-functions at the 
fixed point. Knowledge of the fixed point values of the couplings for arbitrary 
N and the corresponding /3-functions provides the value of v by simple algebraic 
manipulations. 



7 Results and comparison with e-expansions 

In this section, we compare our results for the critical indices at the second-order 
phase transition with those obtained earlier by other authors for large values of 
N. Our numerical results are given in figs. 8 and 9 for the anomalous dimension 
r) and the critical index u, respectively. Different methods have been tried in the 
past to compute these critical exponents for large N |5|, [17]]. Since all methods 
(including ours) are thought to give the most accurate values for large N, a 
comparison of the results is best made for iV — ► oo where rj ~ and v — 1 ~ ■h. 
In table 1 we compare the coefficients Nrj(N) and N(u(N) — 1) as obtained by 
various methods. In our approach the different approximations (and the global 
method) converge very well for these quantities in the large- N limit. Our result 



of the anomalous dimension is given by eq. (|27|) . We have not computed an 
analytical expression for the subleading coefficient of v. This would require a 
systematic expansion of the quantities appearing in the last section in powers of 
1/N. However, this coefficient has been obtained by a numerical fit of the data 
for v{N) to a polynomial in 1/N and we find v — 1 — 

Two papers quote results of a large-iV estimate from a computation in fixed 
dimension d = 3. We should emphasize here that the situation for a gauge theory 
is rather different from the pure 2A r -component scalar theory where exact large- 
N results can be obtained. A loop computation in three dimensions requires the 
presence of an effective infra-red cut-off for all fields. In the pure scalar theory this 
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The limit iV — > oo of 




JV[z/(JV) - 1] 


Our result: 

Previous large- N estimates 17]: 
e-expansion at e = 1 around 

a) d = A-e §: 

b) d = 2 + e, CP^" 1 model, order e 2 17]: 

c) d = 2 + e, CP*" 1 model, order e 4 |17|: 


-0.31 
-2.03 

-9 + 0(e 2 ) 


-1.38 
-4.86 

-48 + 0(e 2 ) 
-2 + 0(e 3 ) 
-6 + 0(e 5 ) 



Table 1 



is provided by the scalar mass term if the phase transition is approached from 



the symmetric phase fll8[ . For the gauge model investigated here the photon 
remains massless in the symmetric phase and produces an infra-red divergence in 
the loop expansion. This problem persists in the large- N approximation. Thus 
perturbative calculations in three dimensions have to introduce some infra-red 
cut-off by hand. We give the estimate of such methods H [I7| in table 1. In our 
approach, the infra-red cut-off is always provided by the averaging scale k and 
no additional assumptions have to be made. We find a substantial discrepancy 
of previous large- N estimates from our result for d = 3. Nevertheless, we expect 



that the relatively crude estimates of ref. [[5], [T7| give rather realistic results for 
d ~ 4 and d ~ 2. For d — 4 the infra-red divergence is only logarithmic and 
the exact form of the infra-red cut-off does not matter. For d = 2 there are 
no transverse photon degrees of freedom and the infra-red problem in the loop 
expansion may be circumvented. 

As an alternative to fixed dimension large- N estimates one may attempt to use 
an extrapolation of results in 4 — e dimensions to e = 1 (e-expansion) || 19 



Results of this method for the critical exponents are also shown in table 1. We 
emphasize that the result in order e differs from our result by factors of about 
30 for both Nt] and N(v — 1)! These rather large differences can directly be 
traced to the threshold effects for the massive fluctuations which are neglected 
in lowest order e-expansion. This can be seen by neglecting the masses in eq. 
(27): The leading 1/N coefficient in the anomalous dimension would then reach 
—9.53 and be of the same order as the e-expansion result —9. For d — > 4 (e — > 0) 
the mass terms are small at the fixed point and the e-expansion becomes reliable. 
For d = 3, however, the mass terms are substantial and the e-expansion fails to 
reproduce the leading 1/N coefficient correctly. This is directly related to the 
absence of a second-order fixed point for moderate values of N within the lowest 
order e-approximation. In contrast to the O(N) symmetric pure scalar theory, 
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N 










1 


-0.134 


-0.170 


0.532 


0.583 


2 


-0.0892 


-0.123 


0.639 


0.568 


3 


-0.0672 


-0.0835 


0.715 


0.664 


5 


-0.0453 


-0.0513 


0.798 


0.777 


7 


-0.0343 


-0.0374 


0.843 


0.832 


10 


-0.0252 


-0.0267 


0.883 


0.877 


30 


-0.00909 


-0.00926 


0.956 


0.956 


100 


-0.00281 


-0.00283 


0.986 


0.986 



Table 2 



where the e-expansion gives very good estimates for the critical exponents |20 
we conclude that the e-expansion gives a very misleading picture of the nature of 
the phase transition for small values of N and sufficiently large A/e 2 . 

Finally, in the limit e 2 — > oo the iV-component abelian Higgs model is thought 

to be described by a pure scalar CP N ~ X model [^T|. We notice that towards the 

infra-red e 2 (k) decreases very fast for large N. The effect of 1/e 2 ^ corresponds 

therefore to a relevant parameter in the appropriately perturbed CP N ^ model. 

Nevertheless, if e 2 K* remains very large at the true second-order fixed point of 

the abelian Higgs model, the critical exponents may be close to the ones of the 

QpN-i mo del. The critical exponents of the CP^ -1 model have been computed 

by an expansion around 2 + e dimensions up to order e 4 [ETJ]. We notice a rather 

poor convergence of the e-series for large N and e = 1. The quoted values of 

the leading 1/N coefficients are also tabulated in table 1 and turn out to be 

much larger than our results. At this stage it is not clear if this is a failure of 

the e-expansion for the CP N ~ X model, or if the 1/N coefficients of the critical 

exponents in the abelian Higgs model differ substantially from the ones of the 
C pN-i model 

Let us next turn to small values of N. The rather small values of the 1/N- 
coefficients of r\ and v — 1 allow a rather smooth extrapolation towards small N 
(cf. figs. 8 and 9). This is in contrast to previous computations. In order to see 
the effects of the approximations we have depicted in figs. 8 and 9 solutions of 
two different approximations to eq. ([14]). For N > 10 their difference is negligible. 
However, for N < 10, effects of higher couplings begin to have a sizeable influence 
on the critical exponents r\ and v. As discussed briefly in sect. 5, we do not expect 
a good convergence of the truncation series for small N. This is also apparent 
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from our results for N <2. The most reliable estimate at present is probably the 
one of the ^-approximation. 

In table 2, we give the critical indices of the second-order phase transition for 
different values of N and for different approximations. We notice that for all 
values of N the anomalous dimension turns out negative. This is in contrast to 
pure O(N) -symmetric scalar theories and is due to the gauge field fluctuations. 
Since for N = 2 the CP 1 model is equivalent to an 5*0(3) model which has 
positive 7] we conclude that for small N our critical exponents differ substantially 
from the CP N ~ X model. It is also clear that the critical exponents deviate from 
the 0(2N) scalar theory which is obtained for e 2 = 0. This is true for all N. 



8 Conclusions 



We have presented a non-perturbative analysis of the scaling behaviour of the 
three-dimensional abelian Higgs model for N complex scalar fields. It is based on 
the investigation of fixed points of non-perturbative flow equations which describe 
the scale dependence of the average action. One fixed point governs the second- 
order phase transition (for the appropriate parameter range) whereas the other 
corresponds to the tricritical point where the second-order transition changes to 
a first-order one. We have found a second-order fixed point for all values of N, 
including the case of the superconductor for N — 1 i. This suggests a second- 
order phase transition for type II superconductors with sufficiently strong scalar 
coupling. Such a picture is consistent with lattice studies i |22[ . We have com- 
puted the critical exponents rj and v in dependence on N. For N = 1 we obtain 
r) between —0.13 and —0.17 and v between 0.53 and 0.58. From the scaling laws 
we infer a between 0.25 and 0.4, j3 between 0.23 and 0.24 and 7 between 1.254 
and 1.135. The phase transition from nematic to smectic-A in liquid crystals is 
thought to be modeled by the same universality class as superconductors. Our 
values for the critical exponents may therefore be compared to the measured val- 
ues for the transition in these liquid crystals. The indices a and 7 agree within 
the experimental uncertainties. For v the experiment distinguishes between par- 
allel and perpendicular directions ||, and direct comparison with our value is 
more difficult. 

We find a very smooth transition between the behaviour for large and small N. In 
8 The case N = 1 is studied in full detail in rcf. fjqj . 

9 Lattice studies cannot distinguish between a second-order and a weakly first-order 
transition. 
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particular, the 1/iV-expansion seems to remain valid down to rather small values 
of N. We were also able to find the tri critical fixed point, even though our global 
method could not be used for iV < 4. We have not attempted here to compute the 
crossover exponents associated to this fixed point. This can be done by methods 
in complete analogy to sect. 6. We expect our results for the tricritical fixed 
point to be less reliable than the ones for the second-order fixed point. This 
is connected to shortcomings of the present truncation that will be discussed 
below. An important message from the present work concerns the validity of 
the e-expansion. In contrast to 0(iV)-symmetric scalar theories we find that the 
e-expansion gives for all N a rather poor description of the phase transition in 
the three-dimensional abelian Higgs model. For large N the e-expansion results 
for the leading 1/N coefficients of rj and v — 1 are off by a factor of 30. For low 
N the e-expansion would suggest a first-order transition for all parameter values, 
in contrast to our findings. 

Let us finally discuss the limits of quantitative validity of our results. This mainly 
concerns the truncation (H) which leads to the non-perturbative flow equations 



(|14D , Qlq) and (p7|). This truncation can be viewed as a lowest order approx- 
imation of a systematic derivative expansion of the most general form of the 
average action Tk- The next order should replace the coefficient Zp,k in eq. (H) by 
a function Zp^{p) and similar for the scalar kinetic terms. This will lead to ad- 
ditional contributions to the flow equations for the potential and the anomalous 
dimension 77^ which involve p-derivatives like Z' 0, [DJ. The most important mod- 
ification, however, is probably the effective replacement of e 2 by a p-dependent 
function. Furthermore, it may become necessary to take the effective momentum 
dependence of e 2 into account. All this can be done by exploiting the exact flow 
equation with less severe truncations. In the present truncation, eq. @, the treat- 
ment of the region around the origin (p = 0) seems not very reliable in the SSB 
regime where the potential minimum occurs for po 7^ 0. In particular, we use at 
the origin effectively e 2 (po) instead of e 2 (0). The resulting errors are particularly 
important for small N and we believe that this is the reason why we did not find 
the tricritical fixed point with the global method for N < 4. Similarly, it seems 
probable that no exact scaling solution exists for the full system of flow equations 
(|14|), (|T5|) and (|D]) for small N p3] . Since the problem appears to be mainly 



related to the behaviour at the origin, a local approximation (in p) seems more 
appropriate at the moment for the treatment of the second-order fixed point. A 
definite computation of the parameter range where the phase transition in super- 
conductors is second-order requires an understanding of the mentioned truncation 
problems and the demonstration of a fully scale invariant solution for Ufc(p). 
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Appendices 



A Functional derivative of the average action 

(2) 

For the derivation of the second functional derivative 1^ , evaluated at (p a = 



^/p6 a i and arbitrary we make the following decomposition |L2[ of the fluctu- 
ations 5<f, 5 A: 



SA U = t u + \/a-^=£ , d»t u = 



6<px = -f=(a + iuj) 

6<p a = —j= (a a + itia) , a > 1. (67) 

Concerning the scalar sector we observe one massive mode a, one massless Gold- 
stone mode uj and 2(N — 1) massless pseudo-Goldstone bosons a a , u a , a > 1. 
Working in the Landau gauge (a = 0) one finds 

J d d x54>rW54> = Jd d x{ £ M [Z F)fe P(-d 2 ) + 2e%, fc p]t" 

+£P(-d 2 )£ + a[Z Vtk P(-D 2 ) sy + U' k + 2pU'l\a 
+u[Z^ k P(-D 2 ) sy + U' k }uo + Wl&y/pZ^t^o 
+2Z^ k uoP(-D 2 ) as a 

N 

+ Y / (a a [Z^ k P(-D 2 ) sy + U' k ]a a 



a=2 

l2\ 



+u a [Z tptk P(-D 2 ) sy + U k ]uj a 
+2Z^ k Cu a P(-D 2 ) as a a )}, (68) 



where 



F[A} sy = \{F[A\ + F[-A]), F[A] as = ±(F[A] - F[-A}). (69) 

Putting A^ = (D 2 = d 2 ) and switching to momentum space, yields 1^ [<p, 0] as 
needed for the derivation of the flow equation ([/]) for the average potential. 



B Threshold functions and constants 

We use the following threshold functions and constants: 

1 roo . 

l d = -k~ d dxx*- l d t \n{P{x)) 

2 Jo 
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Jd 



jd>3 



1 I" 00 rl ~ 

-k 2n ~ d / dxxi^dtPix) 

2 J n 



2 

d-2 



4 



-A- 



4— d 



cixa; 2 — 
ax 



P(x 



■dtP(x) 



'n>l 



/ dxxz-% In (P(x) + ^A; 2 

-^rjk 2n ~ d (°° dxxi~ l d t (P(x)+uk 2 Y n 
2V± Jo ^ ' 



^n,m(^l>^2) 



A .2(n+m)-d / ^1-1 x 



ft 



ft 





(P(x) + u;ifc 2 ) _n (P(x) +uj 2 k 

k 2(n +m -l)- d r dxx i x 
j 

{p{x)+uj l k 2 y n (p(x)+u 2 k'< 
l k 2 ( n +m -i)- d r dxx i x 

2 Jo 



ft 



dP 



dx 
d-2. 



(P(x)+uo 1 k 2 ) n {P{x)+u 2 ti 



4-d 



d 9 d d i / P(a;) — x 

dxx2" 2 - — -In 1 + — ^ - 

cfecft V fc 2 (l + u;) 



(70) 



Here ft stands for acting only on P& contained implicitely in P(x) (see sect. 2). 
These functions, as well as expressions for the asymptotic behaviour of s^(u) for 
uj —>■ oo can be found in [[[]]]. In terms of the above, the threshold function from 
the massive fluctuations for /3 e 2 in eq. ( |PTD reads 

48e 2 K 



s d J2\K,2e 2 K) 



d(d + 2)l d g 
, mft 2 (2A/e,0) 



(d + l)mg 2 (2Ak, 2e 2 /t) - n^ 2 (2A/t, 2e 2 K) 



m d % 2 (0,0) 



(71) 



Furthermore, for large mass terms u 2 » 1 we have 



LJ2 — ►CO 



M^K^) = < o ( Wl ,O)< + 0K (m+1) ) 
lm, l < m K,o; 2 ) = < j0 K,0V 2 - m + O(^ (m+1) ). 



(72) 
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Figure captions 

Fig. 1: The /^-function eq. (|l!J) is given in dependence on A for N = 10 and two 
values of k. The arrows indicate the flow with scale k — > 0. No fixed points 
are obtained when the threshold behaviour is discarded, k = (upper curve, 
right axis). In contrast, the inclusion of the threshold effects with k = 
allows for two fixed points (lower curve, left axis). 

Fig. 2: Eigenvalues of the fluctuation matrix A in dependence on N. The two 
positive eigenvalues indicate that the fixed point corresponds to a second- 
order phase transition for all N. 

Fig. 3: Fixed point value in dependence on N for a ip 4 - approximation. 

Fig. 4: Fixed point value A* in dependence on N for a ^-approximation. 

Fig. 5: Fixed point value e 2 in dependence on N for a (^-approximation. 

Fig. 6: Values for e 2 and A at the second-order fixed point in dependence on N for 
different approximations. 

Fig. 7: Values for e 2 and A obtained with the y? 4 truncation and the global method 
for the tricritical fixed point in dependence on N. The solid line represents 
the large- N estimate given by eq. (|4C|). 

Fig. 8: The critical index rj as a function of N for two different approximations and 
in comparison with previous estimates. 

Fig. 9: The critical index v as a function of iV for two different approximations and 
in comparison with previous estimates. 
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